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Abstract 

In many cases, the critical-state theory can be treated as a sum 
ciently accurate approximation for the modelling of the magnetic prop- 
erties of superconductors. In the present work, the magnetostrictive 
hysteresis is computed for a quite general case of the modified Kim- 
Anderson model. The results obtained reproduce many features of the 
giant magnetostriction (butterfly-shaped curves) reported in the lit- 
erature for measurements made on single-crystal samples of the high- 
temperature superconductor Bi<iSr2C aCu20%. It is shown that addi- 
tion of a contribution to the magnetostriction in the superconducting 
state which is of similar origin as in the normal state, offers a broader 
phenomenological interpretation of the complex magnetostriction hys- 
teresis found in such heavy-fermion compounds as UPts, URu-zSiz or 
UBe 13 . 
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1. Introduction. 

Measurement of the magnetization has been the most commonly used 
experimental technique for the study of flux- pinning. Recently, however, 
magnetostriction studies have been shown informative. For example, Ikuta 
et al. [1] observed giant magnetostriction in Bi 2 Sr 2 CaCu 2 0$ single crys- 
tals in the superconducting state and showed that the largest contribution 
to this effect comes from the interaction between the critical currents in- 
duced by changing the external field and the applied field, i.e., the Lorentz 
force. Much attention has also been directed to dilatometry studies of heavy- 
fermion superconductors. In UPts [2], the irreversible contribution to the 
magnetostriction is small but in UBeis [3], UPd 2 Al 3 [4] and URu 2 Si 2 [5], 
large sample size changes are observed. 

In large fields, when the magnetic induction, B, changes negligibly across 
the sample (e.g., in URu 2 Si 2 , for H > 10 3 Oe [6]), the idealized situation 
for a slab of thickness 2D may be written as B(x) ~ H — (4n/c)j c ■ (D — x), 
with x = at the centre of the slab. Here, H is the external magnetic field 
and j c is the critical current density, (47r/c)j c = —dB/dx. The local Lorentz 
force exerted on the current carriers, / = (l/c)(j c • B(x), is directed towards 
the centre of the sample and leads to its compression, for increasing exter- 
nal field. This force must be compensated by the internal local stress, cr(x); 
da/dx = (1/c) • j c ■ B(x). Hence, with j c = —(c/4tt) ■ B/dx, we find that 
a(x) = — (H 2 — B 2 (x))/8n, for the boundary condition a(D) = 0. For a 
material characterized by the elastic constant en, the relative change of the 
sample size, AD/D, is given by: 



which is equal to —1/(8ttcuD) ■ J (H 2 — B 2 (x))dx. For the present case, 
with (c/An)-j c < B/D, we have: (AD/D ~ -f p /(2-c u )-D. Hence, AD/D 
is directly related to the flux-pinning-force density, f p — (1/c) • B ■ j c . 

2. Magnetostriction in the modified Kim- Anderson model 

For qualitative calculations, we shall consider the simplest possible ge- 
ometry: a slab of thickness 2D, with the external magnetic field H applied 
along its surface. We use a modified Kim-Anderson formula for defining 
the critical current density, j c (B), corresponding to a local gradient of the 
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magnetic induction: 



3c = (b + hy = ~ (c/4?r) ' dB,dx - (2) 

This dependence unifies many forms of the critical state models. The calcu- 
lating the magnetization for this geometry and for j c (B) given by Eq. (2) 
can be found in ref. [7]. This previous work [7] contains an analysis of other 
effects observed in superconductors which can be explained using the model 
of a critical-state. The present work provides additional explanation of these 
phenomena and may be considered as a continuation of our work on this 
subject. 

The calculation of the magnetostriction consists of: i) finding the mag- 
netic field distribution during the magnetization process, ii) integrating the 
equation da/dx = (1/c) -j c -B(x) with the boundary condition that &(D) = 
and Hi) using Eq. (1) for calculation of the magnetostriction. The so- 
lution of the equation (2), determines the flux distribution and has the 
form: B(x) = (±(47r/c) • a ■ (n + 1) • x + /3) 1 ^ n+1 ^ — h, where the inte- 
gration constant f3 is determined from the boundary condition B(±D) = H, 
with the sign of dB(x)/dx dependent on the most recent field change direc- 
tion occurring at point x. For the virgin magnetization curve, f3 is given 
by: p = (H + h) n+1 - H m , where H m = (H* + h) n+1 - h n+ \ and H* 
is the field of the first full flux penetration into the centre of the slab, 
H* = (a ■ (n + 1) • D + Some examples of the flux distribu- 

tion for different situations are drawn schematically in Fig. 1. In this figure, 
Bq is the magnetic induction in the sample centre and is given in Table 1 
for various stages of the magnetization process. The magnetostriction during 
the virgin magnetization process, for H < H*, is found to be; 

87rCll ' AD(iJ) = - ■ ({H + h)^ - h^) 
D H m ^ 

(n + 1) / (H + h) n+3 - h n+3 2h (H + h) n+2 - h n+2 
H m \ n + 3 n + 2 

In order to provide the results for the entire magnetization loop in a trans- 
parent form, we define the following quantities: 
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Figure 1: The magnetic induction distribution in a sample volume when 
a) held increases from to H < H* , b) held increases to H > H* , c) and d) 
when held decreases from H max to H, and e) when held decreases below 0. B 
is the held at the slab centre. B v is the held at the kink of the hux distribution 
at x v , occurring for a certain range of H values after the sweep direction of 
H is changed. x v is equal to 1 - {{H max + h) n+1 - (H + h) n+1 )/{2 ■ H m ). B v 
and B are dehned in Table 1. 

Table 1. Characteristic fields corresponding to the curves a) to e) in 
Fig. 1. By is equal to (0.5 ■ ((Hmax + h) n+1 + (H + - h. 
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Ho,+ 


Hi.- 
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b) 


Bo 


= ((H + h) n+1 -^ m ) 1 /(«+i) -h 
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c) 


Bo 


= ((H max + h)^ 1 - H m ) 1 /(«+D - h 


B v 


B v 
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d) 


Bo 


= ((H + h) n+1 +i/ m ) 1 /(«+i) -h 


H 








e) 


Bo 


= {-(\H\ + h) n+1 + (H* + h) n+1 + / l «+i)V(«+i) 








H 



F iA (H i: _, H h+ ) = h 2 ■ ((H h+ + h) n+1 - (H^_ + h) n+1 ) /H m , 
F h2 (H h _, H it+ ) = -(n+l)/(n+2)-2h-((H h+ + h) n+2 - (H h _ + h) n+2 ) /H„ 
F h3 (H h _, H it+ ) = (n + l)/(n + 3) • ((H h+ + h) n + 3 - (H it _ + h) n ^) /H m , 
with Hi _ and H i + defined in Table 1, for situations corresponding to 
those in Fig. 1. The index i enumerates different branches of the magnetic 
induction curves of continuous j c (x) dependence. Using the above equations, 
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the magnetostriction can be written as 

U i=l,2 

where 

• (Fx if i>+ )'+ F 2 (tf,_, if i>+ ) + F 3 (^,_, if i)+ )) . (5) 

The function S(Hi_, H ij+ ) in Eq. (5) is equal to +1 when Hi - < H ii+ and to 
— 1, otherwise. The remaining part of the magnetostriction hysteresis curve, 
not covered in the Table 1, may be computed using symmetry: AD(—H) = 
Al)(U). 

3. Discussion 

Some examples of the calculated magnetostriction hysteresis curves are 
shown in Fig. 2. A large variety of AD(H) dependences is found: from 
AD(H) proportional to H at H > H\ when n = (Bean model), to AD(H) 
strongly suppressed by field when n > 1. For n > 0, an abrupt passing 
through AD(H) = is found, when the field sweep direction changes. 





Figure 2: The calculated magnetostriction curves, the virgin one and the 
entire hysteresis loop; a) h — 0, while n changes from (the Bean model) to 
n = 2, b) when n = 0.5, h — 1 while H max ranges between 2 and 10. The 
vertical scale is in units of8n ■ c\\ • AD / D. The arrows indicate the direction 
of field change. 
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The field-width of this transition region decreases strongly with an increase 
of the maximum field applied (Fig. 2b). This is due to a low value of 
the critical current density at large fields and a decrease in the field change 
required to reverse the direction of current flow in the sample volume. A 
significant remanent magnetostriction due to a frozen-in magnetic field is 
observed. The features discussed above are commonly observed for high-T c 
and heavy-fermion superconductors. 

An unusual effect has been observed recently [5] for single- crystal of 
URu2Si2- For fields near the upper critical field, H C 2(T), the irreversible con- 
tribution to the magnetostriction changes sign compared to the sign observed 
at lower fields. Moreover, it was observed [5] that the sign of the irreversible 
contribution to the magnetostriction for some experimental configurations 
disagred with the prediction of the critical-state model. On the other hand, a 
large background contribution to the magnetostriction, AD/D = b-B 2 , with 
the parameter b of the order of 10 -16 / G 2 , has been reported for heavy-fermion 
materials [2,3,5], which is apparently of the same origin as the normal-state 
magnetostriction. The role of this contribution on the overall shape of the 
magnetostriction may be considered. We note that the critical-state mag- 
netostriction component results from the calculation of an integral over the 
B 2 (x) dependence. The total magnetostriction in the superconducting state 
is then assumed to be the sum of the critical-state and the normal- state-like 
components. The sign of both contributions must be considered carefuly; 
The critical-state magnetostriction depends on the sign of the magnetic field 
and current while the normal-state contribution depends on the sign of the 
coefficient b. Using these arguments, we may write: 

8tt Cii -AD(H)/D = H 2 - £ (/(ifj _, H i>+ ) - G ■ \ f(Hi_, H i>+ )\) , (6) 

i=l,2 

with f(Hi_,H i>+ ) given by Eq. (5) and G = 8irc n ■ b. Magnetostriction 
calculated from Eq. (6) for a few values of the parameter G is shown in 
Fig. 3. An important feature of these results is that both signs of the 
magnetostriction and both signs of the irreversible contribution to the mag- 
netostriction hysteresis are possible, while the shape of the curves might only 
be weakly affected by the presence of a normal-state contribution (except for 
the presence of a large background proportional to H 2 ). The hysteresis of 
the normal-state-like contribution to the magnetostriction in the supercon- 
ducting state is due to inhomogeneity of the magnetic induction across the 
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sample volume, which is caused by the presence of critical current. It is 
easy to show that at large fields, in particular, near the upper critical field, 
both contributions to the irreversible magnetostriction are proportional to 
the critical current density. Hence, the most intriguing property of the mag- 
netostriction in some heavy-fermion superconductors, i.e. the change of sign 
of the irreversible part of the magnetostriction as a function of the applied 
field, can not be explained by this model, if it is assumed that the parameter 
G in Eq. (6) is independant of the magnetic field value. 
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Figure 3: The calculated magnetostriction curves for n = (the Bean 
model) for several values of the parameter G, determining the normal-state 
contribution to the magnetostriction: G — 2, 1, 0, — 1 and —2, for curves from 
top to bottom, respectively. One must note that the values of G depend on 
the units of the magnetic field assumed, which are here equal to H*. The 
vertical scale is in units of 871 ■ en • AD/D. It has been shifted for some 
curves, for better clarity. The arrows indicate the direction of held change. 

Using Eq. (5), magnetic measurements [6] in conjunction with the value 
of the elastic constant en allow for a determination of AD/D. In URU2S12, 
c n = 2.4 • 10 12 erg/cm 3 [8], yields 2AD/D = 3.6 • 1(T 8 and 2AD/D = 6 • 
10~ 8 , for H\\c and H\\a, respectively, in reasonable agreement with measured 
properties [5]. On the other hand, for UPt 3 , the elastic constant is larger, 
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about 6 ■ 10 12 erg/cm 3 [9], and the critical current density is much lower 
[7], resulting in an irreversible contribution to the magnetostriction which is 
hardly observable experimentally. 

Gloos, et al. [4] reported significant difference between the thermal- 
expansion curves measured during cooling in a magnetic field through the 
superconducting transition and those measured during heating with the same 
magnetic field applied at low temperature, for the heavy-fermion supercon- 
ductor UPd 2 Al 3 . Any hysteresis of this kind finds its phenomeno logical ex- 
planation within the critical state model and is of the same nature as the 
magnetization hysteresis in Mpc(T) and in Mzfc{T) measurements, as de- 
scribed by us in detail in ref. [7]. 

4. Conclusions The present calculations have shown that the total 
magnetostriction in the superconducting state can be separated into critical- 
state and normal-state-like components. In general, this approach has pro- 
vided predictions which are concistent with experimental results in high-T c 
and heavy-fermion superconductors. In particular, the model allows for the 
description of both signs of the irreversible contribution to the magnetostric- 
tion, though does not allow to explain the mechanism causing the change 
of this sign as a function of the magnetic field, which is observed in some 
heavy-fermion systems. 
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